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Abstract 

Neutrinos propagating in matter acquire an effective electromag- 
netic vertex induced by tlieir weak interactions witli the cliarged par- 
ticles in the background. In the presence of an external magnetic 
field the induced vertex affects the flavor transformations of mixed 
neutrinos in a way that, in contrast to the oscillations driven by an 
intrinsic magnetic moment interaction, preserve chirality. We derive 
the evolution equation for this case and discuss some of the physical 
consequences in environments such as a supernova. For small values of 
the square mass difference the resonance for neutrinos and antineutri- 
nos occur within regions which are close. In that case, the resonance 
condition becomes independent of the vacuum parameters and is ap- 
proximately the same for both. 
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As is now well known, under certain favorable conditions, large transfor- 
mations from one neutrino flavor into another may take place in a medium, 
even for small neutrino mixing in vacuum The matter effects on the v 
oscillations are taken into account by means of a potential energy or an in- 
dex of refraction for each neutrino flavor, which can be calculated from the 
background contributions to the neutrino self-energy |Q, |^, ^. In addition to 
the energy-momentum relation of the neutrinos, also their electromagnetic 
properties can be affected in an important manner]^, H, |^, ||, In Ref. 
the general expressions for the electromagnetic form factors of a neutrino in 
an electron gas were calculated and, as an specific application, the correction 
to the index of refraction of a single neutrino in the presence of an external 
magnetic field was determined. In this paper we extend the study of Ref. |0] 
by considering the combined effects of neutrino mixing and the external mag- 
netic field. We derive the equation governing the flavor evolution under such 
conditions and examine the possible effects of strong magnetic fields on the 
neutrino oscillations induced by their effective electromagnetic interaction. 

It is worth stressing the following. The possibility that neutrinos may 
have electromagnetic dipole moment interactions which change left-handed 
neutrinos into right-handed ones can have important consequences in the 
context of the solar- neutrino puzzle 0], and the combined effect of matter 
and magnetic fields on neutrino flavor oscillations and spin precession has 



been studied |]TT1]. However, it is well known that the values of the neutrino 
magnetic moments that are estimated in the Standard Model fall well be- 
low the values for which their effect can be appreciable. While there are 
schemes in which the values of the neutrino magnetic moments lie in the rel- 
evant range, they require ingredients that are not contained in the standard 
Electroweak Theory. On the contrary, the effect we are considering is quite 
different. The neutrino electromagnetic form factors calculated in Ref. |]^ 
preserve chirality and are present even in the Standard Model with massless 
left-handed neutrinos. Those induced terms arise because of the interactions 
of neutrinos with the particles in the background. In the presence of an ex- 
ternal magnetic field the induced form factors, instead of producing spin flip 
transitions, contribute to the index of refraction and modify the resonance 
condition for neutrino oscillations in matter. 

For simplicity we restrict ourselves to the case of mixing between only 
two generations, but the same approach can be applied to more general 
cases. Our starting point is the self-energy of the neutrinos, with momentum 
^ in the presence of a uniform magnetic field. As shown in Ref. it can 
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be written as 

lleff = H + ^ + ($ext)L, (1) 

where is the velocity four- vector of the medium, L = (1 — 75)/2, and 
Bext = l/2e'^'^°'^'Ui^F„/3, where F is the electromagnetic field tensor. In what 
follows we work in the rest frame of the background, where = (1,0) 
k'^ = {u, /?), and i?ext = (0, B) with B being the external magnetic field. 

The coefficients a, b, and c are matrices in the neutrino internal space 
and, in a normal matter gas, composed of electrons, nucleons and their an- 
tiparticles, they are given by 

a = , 

' be 0\ 



b = V2GfQz + 



e 





c = v^G^g'^+(^"j Uj (2) 

where, to order Gp, 

be = V^Gpine-rie) 

with e being the electron charge (e < 0)]] We have introduced the electron 
and positron distributions 

•^t(^) = e/3(B(p)T/^) + 1 ' 
and the corresponding number densities 



He 



(2vr) 



/t(p)- (5) 



In Eq. (1^) we have denoted by Qz and Q'^ the contributions arising from 
the Z-diagram, which are the same for all flavors and, in a normal matter 
background, are irrelevant for oscillations. However, in environments like the 



^We take the opportunity to stress that in Ref. ^ the symbol e stands everywhere for 
the electron charge and not for its magnitude. 
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(b) 



Figure 1: Diagram (a) gives the dominant contribution (of order 
to the neutrino self-energy term of order 5^, arising from the W exchange 
interaction. In the local limit of the W propagator, this diagram is equivalent 
to Diagram (b). 

early universe or the core of a supernova, where the neutrinos represent an 
appreciable fraction of the total density, the neutral-current contributions to 
the potential energy arising from the v-v scattering are not in general pro- 
portional to the unit matrix and should be included in the analysis of the 
resonant flavor transformations |T^. Several approaches exist to describe the 
neutrino oscillations under such conditions, including a Boltzmann-type ki- 
netic approach [|T^] based on a density matrix formalism [H , and a treatment 
based on the FTFT methods|T^. While we are aware that the extra contri- 
butions from the v-u background interactions must in general be taken into 
account in a careful numerical study, in what follows we will not consider 
them further since they can be added at any stage. Of course, they will be 
important in the particular context of the supernova, for example |TB[ . 

In general, the coefficients a, h and c are functions of B. The calculation of 
Ref. corresponds to retaining only the contribution to the coefficients that 
is independent of B. Since we are envisaging a situation where cB could be 
comparable to 6, it is pertinent to ask whether the second and higher order 
terms in B are important or not. In order to answer that, consider the 
contribution to the neutrino self-energy arising from the diagram in which 
two external B lines are attached to the electron line of the M^-exchange loop 
diagram, as shown in Fig. |l]. 

That gives is an additional contribution to such that, instead of the 
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formula for given in Eq. the result is 




(6) 



where 62 is a coefficient independent of B. The crucial point now is to note 
that the term involving the distribution functions f^p cannot be identified 
with the total densities as in Eq. (|[), because the latter quantities must now 
be determined to order also. The diagram for the B^ contribution to the 
electron current density {e'j^e), which can also be written as the trace of 
Sp'j'^ where Sp is the electron propagator, is identical to Diagram (b) but 
with the external neutrino lines removed. Thus, by simple inspection of the 
two diagrams it is easy to recognize that the net number density is given by 
precisely the same factor that appears in Eq. (H), 



This is now the relation that determines the chemical potential in terms of 
the net number density, which is the conserved quantity. If we proceed to 
eliminate the chemical potential in favor of n_ — the result is that the 
formula for be reduces again to that in Eq. (|), while Ce acquires a dependence 
on B through the implicit dependence of the /=p on 5, as a consequence of 
solving for the chemical potential in Eq. (|^). The same argument applies to 
higher order terms in B. This can be seen easily by considering, for example, 
diagrams with more B lines attached to the electron, and noting that they 
are equivalent to the electron loop of Diagram (a) but with the additional 
external B lines. All this amounts to the statement that, to this order in Gp, 
the neutrino self-energy is given by just (e7'^e)7^L. Then, according to Eq. 
([^), be is determined by (e'^e) which is just the net number density. Thus, to 
summarize, the formula for be in Eq. @ is, to order Gp, independent of B. 
On the other hand, by adopting Eq. (j^) as the relation between the chemical 
potential and the net number density, we are neglecting the dependence of 
Ce on the magnetic fieldQ. 

The coefficients b, c in Eq. (0) are written in the flavor basis. In this basis, 
the equation that determines the dispersion relation and wave functions of 

^For the purpose of studying this dependence, a better procedure is probably to start 
from the expression for the neutrino self-energy in terms of the electron propagator in an 
external magnetic field, as recently carried out by Elmfors, Grasso and Raffelt |]l7| 



n 




(7) 
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the propagating modes is 



where, 



m = U 



mi 

777-2 



(8) 



(9) 



with U being the matrix that relates the fields of definite mass to the 
fiavor fields uia = J2i Uaii^Li (« = e, /i). In the Weyl representation 



R 









7] 



(10) 



and Eq. (|^) becomes the set of coupled equations 

{uj — h + a ■ il — ca ■ B)ri — 

{uj — a ■ K,)^ — mi] 



0, 
0. 



(11) 



Using the second of these to eliminate ^ from the first one, yields the following 
equation for r] 



(uj — b + a ■ K — ca ■ B) 



while ^ is then determined as 



m 



7/ = 0, 



= {ui + a ■ k) 



m 

~2 2^- 



(12) 



(13) 



We are interested in the solutions of Eq. ([T^) with positive energy, corre- 
sponding to the (neutrinos) particle solutions. In the absence of the magnetic 
field they correspond to negative helicity spinors of the form 

V = ei,20- , (14) 

where (px is the Pauli spinor with definite helicity, that satisfies 

(CT-«:)0A = A0A (A = ±). (15) 

The ei^2 are vectors in fiavor space, determined by solving the eigenvalue 
problem 

Hei = iViCi , (16) 
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where 



2 

m 



H = K + — + b. (17) 
The expressions for a;i_2 and ei.2 are given exphcitly in Eqs. (A. 10) and (A. 14) 



of Ref. [|15|. To arrive at Eq. (}17\ ), the substitution 



2 2 2 

, ^ m m m 

{u) + a-K)— -7] = — ■ — r]^—r] (18) 

u;^ — uj + K 2k 

has been made. 

In the presence of the magnetic field, the solutions to Eq. (|12|) do not 
correspond to purely negative helicity spinors any longer because the matrices 
a ■ B and a ■ k generally do not commute. Therefore, the solution must be 
sought in the form 

7] = X0_ + , (19) 

where x, x' are two-component vectors in flavor space. Substituting Eq. (|19[) 
into Eq. (|T^) we obtain two coupled equations for x and x', and from them is 
easy to verify that x' ~ {cB ■ kj k)x. Then, retaining terms that are at most 
linear in the small quantities 6, c, rn? /2k in the equation for x, the dispersion 
relations uji^2 and the corresponding vectors x = ei,2 are obtained by solving 



Eq. (16), but with the Hamiltonian 

2 

-> 171 

Hb = K + h- cB ■ k + — (20) 

2k 

in place of H. 

Imitating the arguments given in Ref. , we then arrive at the following 
picture: The Dirac wave function for a relativistic left-handed neutrino with 
momentum k propagating in matter in the presence of an external magnetic 
field is, in the Weyl representation, 

= e^'-' ( ^_ ) X{t) , (21) 

where, in a homogenous medium, x(t) = J2i=i 2(6lx(0))eje~*'^** . For an inho- 
mogenous medium, the flavor-space spinor x(t) is the solution of 

4 = i/.X. (22) 
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which is the extension of the MSW equation to the situation we are consid- 
ering. The components Xe,^i of x(t) give the amphtude to find the neutrino 
in the corresponding state of definite flavor, at a distance r ~ t from the 
production point {to = 0). In Eq. ( pT]) we discarded the right-handed com- 
ponent ipR of the Dirac spinor since it is of order m/2n as compared to the 
left-handed component, and we also neglect the positive-helicity component 
x'(j)^ since it is of order ^ relative to 

We now consider the possible relevance of the effect of the extra terms due 
to the magnetic field on the resonant flavor conversion. At each point along 
the neutrino path, Hb can be diagonalized by the unitary transformation 

U^{r)=( ^^"M , (23) 

with the mixing angle in matter 6m determined by 

Aosin26' , 
sin 2dm = I — , (24) 

J (Ao cos sin 29 

where 

v{r)=K{r)-c,{r)B-k, (25) 

and Aq = {m\ — m^)/2K . The denominator in Eq. (|^) corresponds to 
the difference between the instantaneous energy eigenvalues of the neutrino 
modes u)2 — ^i- For antineutrinos the be and Ce terms change sign and 

- Aosin26' , , 

sin 2^^ = , . (26) 

J (Ao cos 2^^ + ^;)' + A2 sin 29 

According to the above formulas, the mixing angle in matter is modi- 
fied by the neutrino (antineutrino) interaction with the magnetic field. As 
functions of f , sin^ 29m and sin^ 29m exhibit the characteristic form of a Breit- 
Wigner resonance. For the neutrinos the resonance condition (sin2^m = 1) 
is 

f(r^) = Aocos2^, (27) 
while for the antineutrinos we have 

y(r^) = -Aocos2e. (28) 
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The width of the resonance a is given by the length of the interval of values 
of f(r), centered around f(r^), such that sin^ 26^^ > |, and similarly for 
sin^ 26m- In both cases, 

a = 2Aosin20. (29) 



In principle, Eqs. ( |27| ) and ([28| ) can be verified simultaneously in a medium 
with i? 7^ 0. This is a novel feature that contrasts with the situation with- 
out the magnetic field, where only neutrinos, but not antineutrinos, can go 
through a resonant region. For negligible values of Aq cos 26 both conditions 
reduce to 

t;(r^) ~ t;(r^) ~ 0, (30) 

and in such case the resonance for neutrinos and antineutrinos occur within 
regions which are close. More precisely, the above approximations will be 
good to the extent the differences — f| and — f| , where v{f) = 0, are 
small compared with itself. In turn, this implies that 

< , (31) 

where 

1 A, 

(32) 



dv 
dr 



is the characteristic scale length of v at the neutrino resonance, and we have 
taken into account the fact that ~ \\''^r~^r I ■ It is pertinent to remark that 
the condition given by Eq. (|31|) does not necessarily imply that the resonance 
regions of v and V overlap. This will happen whenever the separation between 
the two resonance points is smaller than the average spatial extension of the 
regions. For the neutrinos, the actual spatial extension of the resonance 
region is given by 5y = a/\v'{rR)\ ~ 2A^tan2^, where v' = dv/dr, with 
an analogous expression for 5j7- The condition that both regions overlap is 
< ^i^u+^u), and assuming v'{r^) ~ v'{r^) it simplifies to tan 26* > 1 , 
which requires large values of the mixing angle. 

If the magnetic fields are sufficiently strong that Eq. ( pil| ) can be satisfied, 
then the phenomenon of resonant oscillations can take place, without a severe 
requirement on the masses and mixing angles of the neutrinos. Thus, even if 
the values of these parameters are constrained by the condition for resonant 
oscillations in the Sun and/or other physical phenomena, the supernova may 
simultaneously support resonant oscillations under the conditions just stated. 
In order to estimate the order of magnitude of the magnetic field needed to 
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satisfy Eq. (PPj), in what follows we evaluate Cg in various limiting cases, 
which are easily obtained from Eq. (|]). 

Degenerate gas. In the limit T — 0, we have n-e = and 



-V2GFmeiJ,B ("^) • (33) 



where hb = e/2me is the Bohr magneton. Eq. (|30D then requires 



which can be written in the form 



ussj \0.3j [poj ' ^ ^ 



where p is the mass density, Yg is the fractional number density of 
electrons and p = 10^^ g/cm^. 

Ultrarrelativistic non-degenerate gas. In the ultrarelativistic limit, we 
obtain 

c. = -^^. (36) 

Approximating the Fermi-Dirac distribution by its classical limit we 
obtain 

^ {ue — n-e) = sinh . (37) 
Then, for small values of Ue — n-e, Eq. (|36D reduces to 

Gp 

" - rriefiBP^ine - rie) . (38) 



and the condition in Eq. ( pO]) translates to 



W^gauss J \10MeV 



4x10^) TT^TTITT . (39) 
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Non-relativistic Boltzmann gas. For this case we can borrow the result 
given in Ref. 0, 

Ce = -V2GFf^BP{ne - He) , (40) 
and Eq. ( pUD then becomes 

^ ITf^). (41) 



W^"^ gauss I \WMeV 

It is noteworthy that, in the classical limit, the resonance condition becomes 
approximately independent of the number densities. 

The estimates given in Eqs. (|35|) and (41) indicate that the phenomenon 



we have considered may be relevant in the study of resonant oscillations 
in the supernova, where the densities are typically of order 10^^ g/cm^, and 
magnetic fields of order lO^^gauss have been considered. Our results suggest 
that the effect is worthy of further attention and detailed numerical studies 
in order to asses the range of implications in a concrete way. 

The possibility that the induced neutrino electromagnetic vertex may 
contribute to the index of refraction in the presence of an external magnetic 
field has also been recently pointed out[|18|. However, that work did not 
include the presence of the matter term (be) in the evolution equation, and 
therefore did not consider the implications of the combined effect of that 
term plus the magnetic field term (cg) on the resonant oscillations. 

After this work was completed and while this manuscript was being pre- 
pared, we received a preprint |19| discussing the same effect that we consider 
here. However, the physical outlook and the conclusions reached there re- 
garding the new resonant condition and its implications are different from 
ours. 
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